Giant oscillations of acoustoelectric current in a quantum channel 
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A theory of d.c. electric current induced in a quantum channel by a propagating surface acoustic 
wave (acoustoelectric current) is worked out. The first observation of the acoustoelectric current 
in such a situation was reported by J. M. Shilton et al., Journ. Phys. C (to be published). The 
authors observed a very specific behavior of the acoustoelectric current in a quasi-one-dimensional 
channel defined in a GaAs-AlGaAs heterostructure by a split-gate depletion - giant oscillations as 
a function of the gate voltage. Such a behavior was qualitatively explained by an interplay between 
the energy-momentum conservation law for the electrons in the upper transverse mode with a finite 
temperature splitting of the Fermi level. In the present paper, a more detailed theory is developed, 
and important limiting cases are considered. 



I. INTRODUCTION 

During recent years much attention has been attracted 
by the interaction between surface acoustic waves (SAW) 
and two-dimensional electron gases (2DEG). Most ex- 
periments were performed in GaAs-Al x Gai_ x As het- 
erostructures, and the effects linear in acoustic intensity 

- the sound attenuation and the change of sound velocity 

- were studiedElii These works were aimed to investigate 
the linear response of a 2DEG to a.c. strain and electric 
fields without electric contacts attached to the sample. In 
particular, very interesting experimental studies of these 
effects in-rthe quantum Hall regime were carried out in 
the workam 

The second class of studies deals with the so-called 
acoustoelectric effect, namely a drag of 2D electrons by 
a traveling SAWElij. This effect is due to a transfer of 
momentum from the SAW to the electrons due to SAW- 
electron interaction. As a result, a d.c. acoustoelectric 
current appears in a closed circuit. If the circuit is open, 
the SAW generates a d.c. acoustoelectric voltage across 
the sample. In the simplest case, the current (or the volt- 
age) is proportional to the SAW intensity rather than to 
its amplitude. It is a simple non-linear effect that can be 
also employed to study the 2DEG. 

The first observation of an acoustoelectric current 
through a quantum point contact was reported inE9. 
The authors observed a d.c. current in a quasi-one- 
dimensional channel defined in a GaAs-AlGaAs het- 
erostructure by a split-gate depletion. The conduction 
of such channels is quantized - its dependence on the ap- 
plied gate voltage consists of a set of plateaus divided by 
sharp stepsliinij. 

An important point is that a very specific behavior of 
the acoustoelectric current as a function of the applied 
gate voltage was observed. Namely, the current is not 
quantized. On the contrary, it experiences giant oscilla- 
tions as a function of gate voltage, having maxima near 
the steps between the plateaus of the conductance. 



A semi-quantitative explanation given intJ attributes 
the oscillations to the the electrons in the upper trans- 
verse sub-band of the channel. When the bottom of this 
band is close to the Fermi surface (i.e. near the step of 
conductance) , these electrons at the Fermi level have low 
longitudinal velocity. At some value of the gate voltage 
this velocity can be close to the sound velocity w. Such 
electrons interact strongly with the wave because they 
move synchronously to the latter. At other gate voltages 
there are no electrons which are able to interact effec- 
tively, and the drag is not efficient. 

In the following, we will discuss the acoustoelectric ef- 
fect theoretically in more detail and consider the most 
important parameters and limiting cases. 

II. FORMULATION OF THE PROBLEM 

Consider a ballistic point contact between two regions 
of a 2DEG. According to the Landauer formula, the con- 
ductance of such a junction is determined by transparen- 
cies T n which correspond to different transverse modes 

\nk) = Xn(z,x)exp N J k(£)dA , (1) 

the functions Xn an d the wave-vector k being slowly de- 
pendent on x. Here we label fc the wave vector at x — > oo. 
We have 

*-££|>-(- 2 ^) r .<«->. w 

where /o(enfe) is the Fermi function, while e n k = e n + 
h 2 k 2 /2m is the electron energy for the nth transverse 
mode. Such a concept is based upon the assumption 
that electron thermalization takes place within the re- 
gion ~ £[ n (inelastic relaxation length) near the contact, 
that allows one to reduce the problem to calculation of 
transmittance of non-equilibrium electrons. 
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One can imagine several sources of influence of the 
acoustic wave on the current through the contact. The 
first one is a drag of the 2DEG in the leads. Such a 
drag produces a "phonon wind" inJdic leads which has 
been estimated by Kozub and-Rudir£3 on the basis of the 
model introduced by StredaO. According to this model, 
one can introduce the electron-phonon collision integral 
in the leads, -^^jf (e), as a source of the force acting upon 
the electrons. Then, the current through the contact is 



/ (X 



c-ph 



0. 



£nfc 



T„(e nk ) , (3) 



where (e n fc) is a smooth function, determined by the 
properties of the leads. Consequently, such a contribu- 
tion has the form of a sum of steps-similar to the con- 
ductance that is explicitly stated ir£j. 

However, the piezoelectric interaction between the 
SAW and the electron system in the leads is signifi- 
cantly screened by the 2DEG. Consequently, the steps 
are very small comparing with the effects observed in the 
experiments^, where maxima of the acoustoelectric drag 
current rather than steps were observed. 

As a result, to explain the observed experimental re- 
sults we are left with the region of the quantum channel 
where the piezoelectric field is not too much screened. 
Namely, we analyze the drag due to a momentum trans- 
fer from the SAW to the electron gas inside the quan- 
tum channel. The physical reason of the importance of 
this region in comparison with the much larger region of 
the 2DEG is the absence of effective screening. Such a 
physical picture implies that one has to consider SAW- 
electron interaction inside the channel rather than em- 
ploy Buttiker-Landauer formalism. 



III. GENERAL EXPRESSIONS FOR 
ACOUSTOELECTRIC CURRENT 

In the presence of a harmonic acoustic wave with the 
frequency u> and wave vector q the electrons acquire a 
perturbation 



U 



E 

7ik,n' k' 



k'(lW nk 



O-n'k' 



h.c], 



where 



C n k,n>k>(q) = (nk\exp(iqx)\n'k' 



For simplicity, let us model the quantum point contact as 
a channel with uniform width d. Such an approximation 
is valid if the product of the channel's length L times the 
SAW wave vector q is much greater than 1, qL 3> lJia We 
also assume that the inequality qd <C 1 holds. The last 
assumption allows one to neglect the inter-mode transi- 
tions due to SAW and to take into account only diagonal 
in the mode quantum numbers n, n' contribution. 



In this approximation, C — f3 n (q)5 nn id(k~ k' — q), and 
we can consider the electrons of the n-th mode as moving 
in the effective classical field 

V n (x) = 5R[t//3„ enp(iqx — itot)} . 

Having in mind low enough frequencies, 

hq <C mui,pp , 

one can use the Boltzmann equation for the occupation 
number of the n-th mode, f n k(x) (see, e.g. Ref. |l6|). This 
equation has the form 

d_ y d__ 1 dV n d 
dt dx m dx dv 



+ I)fnk(x)=0, 



where I is the operator describing relaxation of the non- 
equilibrium distribution. We'll specify its form later, tak- 
ing into account two possible relaxation mechanisms: (?) 
impurity scattering, and (ii) escape from the channel. 
Let us expand f n k{x) in powers of the sound amplitude 

U: 



fnk(x) = fo[e nk + V n {x)] + fi + f 2 ■ 



where 



(4) 
(5) 



/i = exp(iqx - iut)] ; fi w <xU 0„ 

The second-order part , f% oc | U (3 n | 2 , is a sum of two items 
- a stationary part, and a part which varies in time with 
the frequency 2uj. We are interested only in the 1st part. 
The second part does not contribute to d.c. current and 
will be omitted. We get from ([!]): 

dfo{e n kY 



Bfi 



oe n k 
dV n (x) gM 
dx dp 1 1 



(6) 



Here we have introduced the operator 

B(q, u) = I + i(qv — uj) 

having the meaning of the operator of the linearized ki- 
netic equation. Angular brackets mean average over time 

t, (• • -)t = {ijJ /2n) Jq 271 ^ ■ ■ ■ dt. Substituting /i and using 
the relationship 

( ^[C u exp(tujt)} exp(iut)] ) t = (l/2)$l(C u D*) 

we obtain 



h = --WP^qu;!- 1 ^ 



{dp 



B~ 



dfo(e n k) 
de n k 



(7) 

As a result, we arrive at the following formal expression 
for the acoustoelectric current (taking into account the 
sum over spins): 



x5R 



n 

dp f d 

2irh \ dp 



dfo{e r 



de n k 



(8) 
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IV. ANALYSIS OF IMPORTANT LIMITING 
CASES 

Eq. (||) is the formal expression for the acoustoelectric 
current. To evaluate it one needs to specify the relax- 
ation operator /, i.e. to discuss sources of relaxation. The 
conventional way to treat relaxation in quantum channel 
is to discuss scattering by individual defects described 
by scattering matrices. However, we want to emphasize 
that the problem of acoustoelectric effect has an impor- 
tant specifics. Indeed, as we will demonstrate later, the 
effect is due to electrons with small momenta k. For these 
selected electrons the effective mean free path, l e g, can 
be much shorter than the mean elastic free path I for 
a bulk 2DEG (see below). We believe that for realistic 
ballistic channels the inequalities £ > L > £ c g hold. As 
a result, one can use a relaxation rate approach, similar 
to the conventional one for bulk systems. The possibility 
for a non-equilibrium electron to escape the channel will 
be allowed for in a model way. 



A. Relaxation rates 

Here we discuss two sources of relaxation, namely an 
elastic scattering within the channel (with the rate 1/r) 
and an escape from the channel due to its finite length L 
(see Fig. 1), an effective rate being \v n \/L. 




this opportunity to emphasize once more that in meso- 
scopic systems an interplay between the impurity and 
other mechanisms of scattering is vety much dependent 
on the problem in question (cf. withE- 1). 

Introducing a correlation function of random impurity 
potential as 

K(r) = (V(r)V(0)) iia 
and its matrix elements 



K mk > tn k = {nk\JC{r)\mk') , 



we get 



T n (k) 



2tt x 

rnk' 



(9) 



If the impurity potential has short range, then \K m k> , n fc| 2 
is independent of the arguments (mk' \nk). Conse- 
quently, 



f 



2n\K\ 2 



-v{e nk ) 



T„(k) h 

where v{e) is the density of states per given spin, 

1 N 1 

n— 1 1 1 

Assuming that the channel behaves as a rectangular box 
with the thickness d in the transverse direction, we get 
N « ^/e F /e*, where e* = h 2 /8md 2 H 

For a wide contact (N — > oo) v = Uf,d, where Vb = 
m/2irh 2 is the 2D density of states. As a result, in a 
channel we arrive at a smooth part of the relaxation rate 
(which is of the order of the momentum relaxation rate 
in the bulk 2DEG, l/rt) plus an oscillating term 
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FIG. 1. Schematic shape of QPC 



Finally, for all the levels except the highest one, one can 
put l/Vb for the impurity relaxation rate, while for the 
highest one the appropriate estimate is 



Here v n is defined as the velocity for the n-th mode at a 
given energy e, 

£n,mv n /h ^ ■ 

For lower transverse modes v n ~ v-p, while for the upper 
one (with the number N) it can be small. On the other 
hand, because of low longitudinal electron velocity, the 
electrons of the upper mode are much more effectively 
scattered by impurities (see below). From this fact we 
come to an important conclusion: for the upper mode 
(which is responsible for the acoustoelectric effect near 
its maxima) the impurity scattering might be important 
even if the contact is ballistic. We would like to use 



7TT b \v N \ 



(10) 



where v* — irh/md. In high-mobility selectively-doped 
structures the scattering potential is smooth. Conse- 
quently, the matrix elements K n k^ n k significantly in- 
crease with the decrease of k. As a result, impurity 
scattering for |£hc upper mode is much stronger than for 
the lower onealj, and the expression (nw acquires an ad- 
ditional large factor a ~ \Kjy : k s[ N,~k~7 K^p ^ | 2 , where 
k s r*> mw/h. 

To take account of the finite length of the contact we 
introduce also the relaxation rate \v n \/ L. As a result, an 
estimate for the relaxation operator of the upper level is 
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max 



1 v* 



Jb\VN\ L 

We observe that there is a border value v c of vn, 



max vf 



a L 
Nil, 



(11) 



(12) 



where both mechanism give contributions of the same or- 
der of magnitude. Here £b is the mean free path in a bulk 
2DEG. At |i>jy| < v c impurity relaxation becomes more 
important than finite size of the contact. 

Expression (|ll]) needs a more detailed discussion. Con- 
sider a ballistic one-dimensional pipe where the parti- 
cles are subjected to a constant force F. As a result, 
they are accelerated as v(t) = vo + (F/m)t, x(t) — 
vot + (F/2m)t 2 . At the time t, the distance between 
a given particle and a particle started after time 9 is 
Sx = 9[vq + (F /m)t\. Consequently, the product of the ve- 
locity at the time t and the local density, v(t) /5x(t) = 1/6 
remains constant keeping constant the current density in- 
side the pipe. 

The total current is determined by the difference be- 
tween the contributions of the particles with opposite di- 
rections of initial velocity. The difference is proportional 
to [F/m)t w FL/mv . Such a result can be reproduced 
by the order of magnitude by the assumption (fill) . 



Here rj = 1/qL has the meaning of the ratio between 
the acoustic wave length and the length of the contact, 
w = uj/q is the sound velocity, while wjvf = vn(^f)- 

We observe that the expression (jlj) consists of the 
product of two functions. The function Ft has sharp 
maxima at — "iVF (see Fig. 2), the width being 



. . k B T 2k B T\ 

d T = mm , \ . 

y mvN V m I 

The properties of F v depend upon rj and the ratio w/v c . 
Here we consider the limiting case 



W/V c < y/fj, 



(15) 



which is relevant to the present experimental situation. 
In this case the impurity scattering dominates for the 
important group of the electrons. If inequality ( |l5| ) is 
met one can put w = 0, and the integrand is symmetric. 
Then 



fix) 



1+X 2 X A + if (I + X 2 ) 2 

Let us consider important limiting cases. 



(16) 



B. Estimates of the acoustoelectric effect 



Now we are ready to make estimates. For simplicity, 
we model the relaxation operator by the following inter- 
polation expression: 



1 V% + V c 

L \vn\ 



(13) 



At vn S> v c it approximately describes escape from the 
channel, while at vp 3> fjv 3> v c it is L-independent and 
equivalent to the inverse life time for the electron with a 
given k due to elastic scattering. We have, 



vnI- 1 = L 



vn\Vn\ 



v 



N 



Integrating Eq. (|^) by parts and taking into account that 
divNl- 1 ) _ 2L \v N \v 2 



dp 



N 



m (v 2 N + v 2 ) 2 



we arrive at the following expression for the oscillating 
part of the acoustoelectric current: 



J 



LeU 2 cu 
irh 



\0 



dv N F v (v N ) F T (v N ) 



F v {v N ) 



rjv 2 N v 2 



(v% + v 2 )[v 2 N (v N - w) 2 + r/ 2 {v 2 N + v 2 ) 2 } ' 
1 1 
Ft{Vn) ~ Ik^T cosh 2 [m(v 2 N -v 2 NF )/ik B T} ' ' ' ' ' 



1. Short waves, rj <Si 1 

At rj <C 1 f(x) has a maximum at x = y/rj with the 
peak value 0.5 and the width ~ y/rj. The shape of the 
oscillations of the acoustoelectric current depends upon 
the relationship between v C y/rj and St- At v C y/rj <C St 
one can replace 



As a result, 



LeU 2 LUy/rj 2 



s/2hv c 



\P N VF T {v c y^). 



(17) 



In this case, we have a peak near wjvf = (i.e. exactly at 
the step), its shape being determined by the derivative 
of the Fermi function, Ft ■ 

At v C y/fj ^> St the function Ft behaves as 

S[(m/2)(v 2 N - v 2 NF )} = (l/mv NF )5{v N - v NF ) . 



Consequently, 



2LeU 2 oj 
irfrvNFmv 2 



\(3 N \ 2 f 



(18) 
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Note that in both cases the current is independent of the 
channel's length L. It is natural, because we consider the 
situation where the intra-channel impurity scattering is 
the most important relaxation mechanism. 

At w/v c 3> ^/rj the escape of non-equilibrium electrons 
is most important. In this limit one obtains j oc L.Eil 
Consequently, measurements of the L-dependence of the 
acoustoelectric current would help to discriminate be- 
tween different relaxation mechanisms. 



2. Long waves, rj 3> 1 

If the wave length of the SAW is greater than the chan- 
nel length, the function f(x) can be approximated as 



f(x) = rf 



(l + .x 2 ) 



2^3 



(19) 



It has a maximum at x = 1/^/2 and width of the order 
1. Consequently, in dimensional variables the width is of 
the order of v c . It has to be compared to the difference 



S v = \v N - v N -i\ w ^4Ne*/m w v F N~ 1/2 



If 



> S v 



or L > I 



the oscillations are not pronounced. On the contrary, 
at L <C I the oscillating part is pronounced. Again, 
its shape is determined by the relationship between the 
widths of the functions F v and Ft- At v c St the 
result is given by Eq. ( |l8|) with f(x) taken from (19). 
In the opposite limiting case, f(x) can be replaced by 
(n/16r])S(x — 1/^/2). Consequently, the shape is deter- 
mined by the function Ft, like in Eq. ( p"7| ) 



• The contact has finite length L, the corresponding 
rate being ~ \v\/L. This rate decreases as v — > 0. 

• Density of states in a QPC is an oscillating function 
of the energy, it diverges at the thresholds corre- 
sponding to the filling of new levels. Consequently, 
the disorder-induced relaxation rate for the upper 
mode increases at the threshold. 

As a result, the total relaxation rate is a non-monotonous 
function of the electron velocity. The analysis given 
above leads to the conclusion that coupling is optimal 
for the electrons having their velocities in a relatively 
narrow range, the central velocity s ~ v c ^/rj being small 
comparing to the Fermi velocity . Above we have esti- 
mated the position of the center and an effective width 
of the important region which is shown in Fig. 2. 




Velocity 



FIG. 2. Scheme of energy and momentum conservation 
laws for electron-SAW interaction 



V. DISCUSSION 

Let us discuss qualitatively the picture of the acousto- 
electric effect. The linear response of the electrons to the 
SAW with a given wave vector q is proportional to the 
effective "interaction time" (qv — uj)~ 1 , during which an 
electron with the velocity v moves in an almost constant 
field. At small v, or near the resonance (v = w), this time 
diverges, and relaxation becomes important. In fact, the 
coupling is proportional to ^ qv _^ 2+I2 ■ Consequently, to 
get an effective coupling both the electron velocity and 
the scattering rate have to be small. 

In a homogeneous 2D system, the most important re- 
laxation mechanism is disorder-induced scattering, the 
scattering rate being proportional to the 2D density of 
states. Thus, the electron-SAW coupling for such a sys- 
tem is determined by the product qt. In a point contact 
relaxation differs from the case of homogeneous 2DEG 
due to following reasons. 



On the other hand, electronic states can contribute to 
interaction only if their energies are in the vicinity of the 
Fermi level. Consequently, for a given mode number n 
the electron velocities have to belong to a narrow interval 
centered around some velocity v/vf- One can move these 
intervals by changing the gate voltage. 

To get a non- vanishing contribution to the current, 
one of these regions has to overlap with the region cen- 
tered around s (see the Fig. 2), which is possible only for 
the upper level. As a result, the acoustoelectric current 
experieaces giant oscillations as a function of the gate 
voltage.c3 Indeed, in the plateau region for the conduc- 
tance the velocity vnf is large enough and cannot overlap 
with the hatched region near s. Then, as the system is 
driven to the step, vnf decreases and important regions 
start to overlap. Consequently, the current increases. 

The fine structure of the peaks needs more careful dis- 
cussion. The current theory can lead to quantitative re- 
sults only when impurity relaxation is more important 
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than the escape from the channel. In the opposite lim- 
iting case the model assumption (^) can provide only 
order-of-magnitude estimates. The exact results in such 
a situation depend both on the length and the shape of 
the channelEil. 

An important source of the fine structure of the peaks 
could be a gate-voltage-dependent screening of the cou- 
pling. This problems needs much more careful considera- 
tions. For a rough estimate, we can introduce a screening 
factor as 



a rough estimate can be given 



U 2 = 



u 2 



i + g[v(q,u)/s\ 



where a(q, to) is the effective conductance while g is a 
(small) geometry-dependent dimensionless factor. The 
effective conductance can be estimated as a sum of 
the contributions of lower transverse modes (which is a 
smooth function of the gate voltage) and of the contri- 
bution of the upper mode 



U ~ yJxSe 2 /uj. 

Here x is piezoelectric coupling constant which can be de- 
termined by experiments on the change of sound velocity 
due to 2DEG. To check the presented theory quantita- 
tively it seems important to measure acoustic intensity 
independently. 
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The latter can be estimated near the peak of the drag 
current in the same way as for a long wire because usu- 
ally L/s > t. Considering a piece of wire with the length 
w q^ 1 , we get the estimate C ~ q^ 1 for its capacitance 
(up to logarithmic terms), and R ~ qai(q, oj) for its resis- 
tance. Here <j\ is the ID conductance calculated from the 
kinetic equation. As a result, the dimensionless screen- 
ing parameter can be estimated as (loRC)^ 1 ~ qo\/s. 
Consequently, the screening parameter is 



cr osc (q,uj)/s = y-^ I dp [ vB ' 



1 



dfo{e n k) 
de n k 



sh qv n T 



(20) 



Far from the maximum of the acoustoelectric current 
this quantity appears small. However, near the maxi- 
mum it can be of the order 1, leading to a decrease of 
the current. Probably it is the origin of the double-peak 
structure o£4he 1st peak of the acoustoelectric current, 
observed ir£3. With the increase of the number of occu- 
pied modes the smooth part of screening increases. As a 
result, the current oscillations' amplitude must decrease, 
and the double-peak structure has to be less pronounced. 
Such a behavior is in a qualitative agreement with the 
experiments. 

Unfortunately, it is very difficult to give realistic esti- 
mates for the coupling constant connecting the intensity 
S of SAW and the amplitude U of the electron's Fpten- 
tial energy. According to the experimental rcsultsEj we 
believe that it is determined by piezoelectric interaction 
iry-the channel. Otherwise, under the conditions realized 
irJ13 steps predicted inE3 would be observed rather than 
giant oscillations. However, the SAW in a layered struc- 
ture has a complicated polarization structure, and only 
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